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■ 1 Introduction 

' In the Fault- Tolerant Facility Placement problem (FTFP), we are given a set T of sites at which 

facilities can be built, and a set C of clients with some demands that need to be satisfied by different 
facilities. A client j has demand rj. Building one facility at a site i incurs a cost /j, and connecting 
one unit of demand from client j to a facility at site i £ J costs dij. The distances dij form a 
metric, that is, they are symmetric and satisfy the triangle inequality. In a feasible solution, some 
number of facilities, possibly zero, are opened at each site i, and demands from each client are 
connected to those open facilities, with the constraint that demands from the same client have to 
be connected to different facilities. Note that facilities at the same site are considered different. 
It is easy to see that if all rj = 1, then FTFP reduces to the classic uncapacitated facility 
[~^. | location problem (UFL). If we add a constraint that each site can have at most one facility built, 

then we get the Fault-Tolerant Facility Location problem (FTFL). One implication of the one site 
per facility restriction in FTFL is that maxj g c r j < \J~\, while in FTFP rj can be much bigger than 

m 

UFL has a long history and there has been great progress in designing better approximation 
algorithms in the past two decades. Since the publication of the first constant approximation 
algorithm by Shmoys, Tardos and Aardal [10], we have seen a number of techniques that are 
applicable in devising approximation algorithms for UFL with good approximation ratios. Using 
LP-rounding, Shmoys, Tardos and Aardal [10] obtained a ratio of 3.16, which was then improved by 
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Chudak [4J to 1.736, and later by Sviridenko [TT] to 1.582. Byrka [2] gave an improved LP-rounding 
algorithm, combined with the dual-fitting algorithm in [6], achieving a ratio of 1.5. Recently Li [9] 
showed that with a more refined analysis and randomizing the scaling parameter, the ratio can be 
improved to 1.488. This is the best known approximation result for UFL. Other techniques include 
Jain and Vazirani's [7] primal-dual algorithm with ratio 3, and Jain, Mahdian, Markakis, Saberi 
and Vazirani's |6j dual-fitting algorithm with ratio 1.61. Arya et al. [1] showed that a local search 
heuristic achieves a ratio of 3. 

FTFL was first introduced by Jain and Vazirani [8] and they adapted their primal-dual algorithm 
for UFL to obtain a ratio of 31n(maxj g c r j)i which is logarithmic in the maximum demand. Guha, 
Meyerson and Monagala [5J adapted the Shmoys et al. 's [10J algorithm for UFL to obtain the first 
constant approximation algorithm. Swamy and Shmoys [12] improved the ratio to 2.076 using the 
idea of pipage rounding. Most recently, Byrka, Srinivasan and Swamy [3] improved the ratio to 
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1.7245 using dependent rounding and laminar clustering. All the constant approximation algorithms 
are based on LP-rounding. 

FTFP was first introduced by Xu and Shen [13] as a generalization of UFL. They extended 
the dual-fitting algorithm [6] to give an approximation algorithm with a ratio claimed to be 1.861. 
However the algorithm runs in polynomial time only if maxj^c r j is polynomial in 0(|.F| • |C|), and 
the analysis of the performance guarantee seems flawed. 

Our approach is to reduce FTFP to FTFL. We then use Byrka et al. 's [3] result on FTFL to 
achieve the ratio 1.7245. 



2 The LP Formulation 

The FTFP problem has a natural IP formulation. Let yi represent the number of facilities built at 
site i and let represent the number of connections from client j to facilities at site i. If we relax 
the integral constraints, we obtain the following LP: 



minimize J2ieF fm + Eie^jee d v x v i 1 ) 
subject to yi — > 0, Vi € J 7 , j € C 

Xij >0,yi> 0, VieTJeC 



The dual program is: 

maximize Yljec r J a j (^) 
subject to Ylj&cPij ^ fh Vi G J 7 

a-j - Pij < dij, Vi G T,j £C 

aj>0,/%>0, yi£T,j£C 



3 A 1.7245- Approximation Algorithm 

In this section we give an algorithm based on a reduction to FTFL. Moreover, we show that if 
we have a p- approximation algorithm for FTFL, then we can use it to build an approximation 
algorithm for FTFP with the same ratio. To be more precise, we assume that p is the ratio of the 
cost of the FTFL solution to the fractional optimal solution cost of the natural LP formulation for 
FTFL. 

A naive idea is to split the sites so that each site is split into P = maxj g c rj identical sites. 
Now we can restrict each split site to have at most one facility. Since we never need more than P 
facilities open at the same site in the original instance, the FTFL instance after split is equivalent 
to the original FTFP instance. However, since P might be large, this reduction might result in an 
instance with exponential size. 

1 Adding a constraint y; < 1, Vi G T to the LP (TI]) results in the LP for FTFL. 
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However, we can modify the naive approach to obtain an FTFL instance with polynomial size, 
using an optimal fractional solution to the LP ([I]). 

Let (x*,y*) be the fractional optimal solution to LP ([1]). Also let 

jji = max{0, [y*\ - 1}, = min{ [x*j\ , 

and 

= x^j X{j, Di = Di Vi- 

Let X be the original FTFP instance. We define an FTFP instance X\ with the same set T of sites 
and C of clients, the same distances dij and the same facility costs /j, except that the demands are 
fj = Y^ieF^ij- Let &i be the solution defined by (xij,y~i). Clearly Si is a feasible integral solution 
to the instance T\. 

Another FTFP instance X2 is defined using the same parameters except that the demands are 
fj = rj — fj = Yi^jrXij. We claim that (xij,jji) is a feasible fractional solution to the instance 
I2 El- To ensure feasibility, we only need to show Xij < fji for all i € T, j € C. The proof proceeds 
by considering two cases. Case 1 is that y* < 1. Then we have x^ = x*j < y* = fji as needed. Case 
2 is y* > 1. From the definition we have fji > 1. We further consider two subcases: if x^ = fji, 
then x^ = x*j — x^ < y\ — jji <yi and we are done. The other subcase is x^ < fji. It follows that 
< Xij < 1. Again we have x^ < 1 <fji- The above takes care of all cases, so we have shown the 
feasibility of the fractional solution (xij,y~i) to the instance Xi- 

One more observation about the instance X% is that fj < 2n where n = |.F|, since fj = YlieF^ij 
and for every i S J, j G C, we have < x^ < 2. Therefore we have P = max Jg cfj < 2n. Now 
that we have the demands being polynomial w.r.t. the input size of the original FTFP instance, we 
can split each site into P new sites and treat the instance as an FTFL instance. Solving the FTFL 
instance with a ^-approximation algorithm, we have an integral solution S2 with cost no more than 
p ■ LP*(Z2), where LP*(Z2) is the cost of a fractional optimal solution to the instance X2. 

Combining S\ and 52, we have a feasible solution to the original FTFP instance. We now argue 
that the solution is within a factor of p from the cost of an optimal fractional solution to the LP ([I]) . 

First we observe that cost(Si) = Yli^F fiVi + YlieJ 7 jecdij^ij- Secondly, since (x~ij,yi) is a 
feasible fractional solution to the instance I2, we have X^eJ 7 fiVi + Sie-FjeC dijXij > LP*(Z2). Our 
solution has total cost 

cost(Si) + cost(S 2 ) < ^2 hm + ^2 dijXij + p • LP*(X 2 ) 

< ^2 + ^2 di i^ + p ^2 foy* + ^2 dijxij) 

^ ~] fiVi ~l~ ^ ] dij Xij) 

= p(^M + %4) = p-LP*(J) 

< p • OPT(X) 



Note that if we take Xij — \x*j\ , fji = [y*\ and define x lj ,y i , r 3 in a similar way, then {xij,yi) may not be feasible 
to the instance with demands {fj} because we might have Xij > tji if < x*j < 1 and y* = 1. 
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The currently best known approximation algorithm for FTFL achieves a ratio of 1.7245 [3], so 
our algorithm is a 1.7245-approximation algorithm for FTFP. 

4 Approximation for Large Demands 

If we seek an approximation algorithm with ratio as a function of R = mmj^c r j, then we show 
that we can achieve a ratio of 1 + 0(n/R) where n = \J-\. The hidden constant in the big-0 is the 
same as the best approximation ratio for FTFL, which is known to be at most 1.7245 [3]. H 

Our algorithm uses the 1.7245-approximation algorithm for FTFL as a subroutine. First we 
solve the LP (JTJ) and obtain (x*j,y*) as the optimal fractional solution. Now we round down the 
fractional values as before and the rounded down values Xij = [x*j\ , y-i = [y*\ form part of the 
final solution, denoted as Si. Now each client j has a reduced demand fj = rj — Y^ieT^ij — n ~ 1) 
where n = Therefore we can solve this instance with reduced demands, denoted as Z 2 , by 
solving an equivalent FTFL instance: We replicate each site with n — 1 copies since no site needs 
to have more than n — 1 facilities open. Now we can restrict each duplicated site open at most 
one facility and this gives us an FTFL instance with polynomial size to work with. We then use 
the 1.7245-approximation algorithm to solve this FTFL instance to obtain an integral solution S 2 
with cost(S 2 ) < p ■ LP*(Z 2 ) with p = 1.7245. We now show that LP*(Z 2 ) < ^LP*(Z) < ^OPT(X), 
where X is the input FTFP instance. To see this, consider three FTFP instances: X\ is the instance 
with all rj = R, and Z 2 is the instance with all rj = fj, and I3 is the instance with all demands 
rj = n — 1. We have LP*(Z3)/LP*(Zi) = (n — 1)/R since any fractional solution with cost Z of an 
FTFP instance with uniform demands s can be scaled down by s to obtain a feasible solution for 
the same instance with demands set to 1 and the cost being Z/s. We also have LP*(T3) > LP*(l2) 
since every demand in Z3 is at least as large as that in Z 2 . Similarly LP*(Z) > LP*(Xi). Therefore 
we have 

LP*(X 2 ) LP*(Z 3 ) = n-1 n 
LP*(X) - LP*(Zi) R ~ R' 

By combining S\ and S2 we have a feasible solution to the original FTFP instance. Now we 
bound the total cost. First notice that cost(S\) < LP*(X) < OPT(Z). By the previous argument, 
we have cosi(5 2 ) < p^OPT(I). Therefore, the total cost of our solution is 

Tl 

cost{S!) + cost{S 2 ) < (1 + p-)OPT(X), 

R 

which is (1 + 0(^))OPT(I). The approximation ratio approaches 1 when R becomes much larger 
than n. This is perhaps not too surprising as when all the demands are large, the fractional optimal 
solution to the LP ([I]) is close to integral relative to the demands, so their costs are close. 
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